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Abstract. We discuss high-resolution, time-resolved spectra of the caustic exit of the binary microlensing event OGLE 2002- 
BLG-069 obtained with UVES on the VLT. The source star is a G5III giant in the Galactic Bulge. During such events, the source 
star is highly magnified, and a strong differential magnification around the caustic resolves its surface. Using an appropriate 
model stellar atmosphere generated by the PHOENIX v2.6 code we obtain a model light curve for the caustic exit and compare 
it with a dense set of photometric observations obtained by the PLANET microlensing follow up network. We further compare 
predicted variations in the Ha equivalent width with those measured from our spectra. While the model and observations 
agree in the gross features, there are discrepancies suggesting shortcomings in the model, particularly for the Ha line core, 
where we have detected amplified emission from the stellar chromosphere after the source star's trailing limb exited the caustic. 
This achievement became possible by the provision of the very efficient OGLE-III Early Warning System, a network of small 
telescopes capable of nearly-continuous round-the-clock photometric monitoring, on-line data reduction, daily near-real-time 
modelling in order to predict caustic crossing parameters, and a fast and efficient response of a 8m class telescope to a "Target- 
Of- Opportunity" observation request. 

Key words. Techniques: gravitational microlensing - Techniques: high resolution spectra - Techniques: high angular resolu- 
tion - Stars: atmosphere models - Stars: individual (OGLE 2002-BLG-069) 



1. Introduction 

Near extended caustics produced by binary (or multiple) lenses, 
the source star undergoes a large total magnification in bright- 
ness. Furthermore, its surface is differentially magnified be- 
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* Based on observations made at ESO, 69.D-0261(A), 269.D- 
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cause of a strong gradient in magnification. The relative lens- 
source proper motion is typically slow enough to allow the light 
curve to be frequently sampled. This translates to a high spa- 
tial resolution on the source star's surface and hence permits its 
radial brightness profile to be inferred from the observations. 
Over the past four years, coefficients characterizing linear or 
square-root limb-darkening profiles have been obtained with 
the microlensing technique for several Bulge giants and sub- 
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giants jFields et alj2003l and references therein) and a main se- 
quence star lAbe et al.l2003b. A new generation of stellar atmo- 
sphere models JOrosz & HauschildttZOOOb have revealed limb- 
darkening laws that a re significantly different from the tra- 
dition al analytic ones dBrvce et alJl2002l Iciaret & HauschildJ 
120031) . The centre to limb variation of spectral lines can show 
markedly different behaviour from that of the continuum. Most 
moderately strong and weak lines weaken towards the limb, but 
reso nant scattering lines can vary in a much more pronounced 
way JLoeb & Sasselovll995b . In cool giants, H/3, being formed 
lower in the atmosphere, is more limb-darkened than Ha. TiO 
and Can show strong variations, and some lines and bands 
may even be limb-brightened. Intensive spectroscopic moni- 
toring of a caustic exit at high resolution with high S/N should 
reveal temporal variations in the equivalent widths of promis- 
ing spectral lines that can be compared with pred ictions from 
stellar atmosphere models jHevrovskv et all2000l) . 

We present here the first photometric and spectroscopic 
monitoring campaign that has successfully been performed at 
high r esolution with dense sampling. Previously. ICastro et alJ 
fcOO lb obtained two KECK HIRES spectra of EROS 2000- 
BLG-5, but missed the trailing limb of the caustic where the 
effects are str onger, while the data in the qualitative analysis of 
lAlbrow et al.1 (2001) involved dense coverage but at low resolu- 
tion. lAfonso et al.ld200lb provided a model that reproduces the 
photometric data and found an excess of Ha at the limb. They 
attribute the excess to chromospheric emission, but it could be 
due to shortcomings of the synthetic spectra. Using a determi- 
nation of the spectral type of the source star, we have com- 
puted the limb-darkening profiles in R and I from appropri- 
ate PHOENIX synthetic spectra of the source star and fitted 
a fold-caustic model to our photometric data obtained during 
the caustic exit. This model has been used here to compute 
the synthetic spectrum for wavelengths around the Ha-line. A 
more detailed report on the determination of the stellar param- 
eters and t he analysis of other s pectral lines will be presented 
elsewhere JBeaulieu et all2004T) . as will the details of a full bi- 
nary lens model fi t to the PLANET observations of the event 
jKubas et alJ2004 . 

2. OGLE 2002-BLG-069 photometry and 
spectroscopy 

The PLANET collaboration, comprising six different tele- 
scopes, namely, SAAO 1.0m (South Africa), Danish 1.54m 
(Chile), ESO 2.2m (Chile), Canopus 1.0m (Australia), Mt. 
Stromlo 50in (Australia) and Perth 0.6m (Australia), com- 
menced photometric observations of OGLE 2002-BLG-069 
alerted by the OGLE-III early Warning System dUdalskil2003b 
in early June, 2002. From online data reductions on 25 June 
it was realised that the event involved a binary lens and a 
bright giant source star. With the source taking ~ 1.4 days 
to cross the caustic on entry, it appeared to be an excellent 
candidate for time-resolved spectroscopy of the caustic exit. 
Using the predictions based on modeling our photometry, and 
thanks to excellent coordination with the staff at La Silla and 
Paranal, very good coverage was obtained during the caustic 
exit together with post-caustic reference spectra. These obser- 



vations were performed using the U VES spectrograph mounted 
on Kueyen (VLT UT2) as part of Target-of-Opportunity and 
Director's Discretionary time. Thirty-nine spectra with expo- 
sure time of 20 min (1 hour for the post caustic observa- 
tions) were obtained alternately in the so-called standard set- 
tings, 580 and 860, covering the full 4800 - 10600 A range 
at a resolution of 30000 with S/N ranging from 50 to 130. 
From an analysi s of the curve-of-gro wth of 100 Fe I and Fe II 
lines as done bv lMinniti etaD (12002b . and independent ly from 
the study of the Call and Mg lines diOr gensen et alJll992b . 
a good fit to all the data was obtained with a plane-parallel 
model atmosphere having T c g = 5000 K, log(g) = 2.5, 
fWb = 1-5 kms -1 and [Fe/H] = —0.6, The calibration of 
the MK spectral type gives, for a G5III star, T cS = 5050 K, 
M v = +0.9, (V - I) = 0.95, M/Mq = 1.1 and R/R Q = 10. 
At HJD=245 2455. 2754, the source was amplified by 17.7. The 
measured color from SAAO data are I = 13.01 ± 0.01, 
(V-I) = 2.06±0.02,soE(V-I) = 1.11 ±0.02. We adopt a 
conservative Ay /E(V — I) = 2.2 ± 0.3, and then derive a dis- 
tance of 9.4 ±1.4 kpc. 

3. Modeling of the caustic exit 

With p denoting the fractional radius of the source star, 
we approximate its wavelength-dependent brightness profile 
I x (p) = I£, x (p) by concentric rings of constant intensity. A 
model atmosphere for the stellar parameters given in the previ- 
ous section was computed with the PHOENIX grid v2.6 (based 
on the code described bv lAllard et alJ J200lb but using spher- 
ical geometry and spherically symmetric radiative transfer). A 
synthetic spectrum was calculated from this model at a spectral 
resolution of 0.05 A at 128 steps in source radius. The inten- 
sity profile was interpolated with cubic splines at 1000 points 
equally spaced in cos 1? = y(T~— p 2 ) where ■& is the emergent 
angle, so that the width of the rings of constant intensity de- 
creases toward the stellar limb. For a broadband filter (s), the 
intensity profile ^ s > (p) is obtained by convolving the source 
brightness profile with the filter, CCD transmission and atmo- 
sphere response functions. The observed magnitude m^ s '(t) 
reads m^(t) = m { s s) - 2.5 lg {A^ (t) + g^} where m { s s) 
is the intrinsic source magnitude, g^ — Fg /-Fg is the ratio 

(s) (s) 

between background flux and intrinsic source flux Fif>, 
and (t) the source magnification at time t. 

The photometric data presented in the lower panel of Fig.^ 
show a peak with the characteristic shape resulting from a fold- 
caustic exit. For a source in the vicinity of a fol d caustic, the 
magnification A^ (t) can be decomposed (e.g., lAlbrow et alJ 
Il999ab as A^(t) = A ( c f it (t) + A othcr (t), where A^ t (t) de- 
notes the magnification of the critical images associated with 
the caustic, and A ot hcr(t) denotes the magnification of the re- 
maining images of the source under the action of the binary 
lens. 

Let us consider a uniformly moving source crossing the 
caustic within 2 At and exiting the caustic by its trailing 
limb at tf. If one neglects the curvature of the caustic and 
the variation of its strength over the size of the source, 
the magnification of the critical images reads A^) lt (t) = 
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Fig. 1. Equivalent-width variation of the Ha core (top panel) 
and wings (second panel from the top); the open circles in the 
plots of the equivalent-width variation correspond to our U VES 
data, while the adjoining solid lines represent the model pre- 
dictions over the course of the caustic passage. The big dots 
correspond to the spectra of 10 July, UT 02H58. Third panel: 
model light curves and photometric data, where R-band data 
are plotted in red and I-band data in blue; model residuals from 
the chosen PHOENIX atmosphere (fourth panel) and from the 
linear limb-darkening law (bottom panel), with the same color 
convention. The model parameters can be found in Tab.Q The 
majority of model residuals are below the 2 % level. The post 
caustic observations have been plotted at 2465.85, but have 
been taken on 16 August. 

q crit Gf ( — -^l 1 ; where Gf(n\ is a characteristic 

function JSchneider & Wagoner! 1 1987D which solely depends 
on the intensity profile £w. For \uj(t — t{)\ -C 1, the mag- 
nification of the non-critical images can be approximated by 
ArtherC*) - aothcr [1 + u(t - tf)] where a ot her is the magnifi- 
cation at the caustic exit at time tf and u> measures its variation. 

As the period during which our generic fold-caustic model 
is believed to be a fair approximation, we choose the range 
2463.45 < HJD - 2450000 < 2467. Restricting our atten- 
tion to those data sets with more than two points in this region 
leaves us with 29 points from SAAO and 15 points from UTas 
in I as well as 98 points from the ESO 2.2m in R, amounting to 
a total of 142 d ata points. From a b inary lens model of the com- 
plete data sets jKubas et al.l l2004) for these observatories and 
filters, we have determined the source magnitudes ml and 
the blend ratios g^ s \ With the adopted synthetic spectra, we 
compute the stellar intensity profiles for I and R, and use these 
for obtaining a fit of the generic fold-caustic model to the data 
in the caustic-crossing region by means of x 2 -minimization, 
which determined the 5 model parameters tf , At, a c rit> Mother, 
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Table 1. Model parameters of a fit to data obtained by PLANET 
with the SAAO 1.0m, UTas 1.0m, and ESO 2.2m. on the mi- 
crolensing event OGLE 2002-BLG-069. While source magni- 
tudes and blend ratios </M result from a fit involving a bi- 
nary lens model to the complete data sets, the remaining param- 
eters have been determined by applying a generic fold-caustic 
model with PHOENIX limb darkening solely to the data taken 
around the caustic exit. 



Parameter 


Value 


Parameter 


Value 


mf A0[7] 


16.05 


tf (days) 


2465.637 




16.05 


At (days) 


0.7297 


ESO 2.2m [Hi 


16.3 




19.60 


ff SAAO [I] 


0.16 


Mother 


7.011 


ff UTas [/] 


0.064 


u> (days) 


-0.04519 


^ESO 2.2m [R] 


0.0083 







and w as shown in Tab.^ If a classical linear li mb darkening 
is inclu ded in the list of parameters to fit (as in lAlbrow et alJ 
Il999bh . the best fit is obtained with T = 0.5. 

The modeled light curve and the structure of the residu- 
als of the two fits are given in the lower panels of Fig.^ The 
fit with adopted PHOENIX limb darkenings clearly show sys- 
tematic trends at the level of 1-2 %. The same trends in the 
residuals can be seen from the figure, which suggests some spe- 
cific features lying in the real stellar atmosphere. We checked 
whether this discrepancy can be caused by the straight fold 
caustic approximation, but from running the already mentioned 
global binary lens model neither the very small curvature of the 
caustic nor the source trajectory can explain this systematic ef- 
fect. 

4. Comparison between synthetic spectra and 
spectroscopic data 

The fold-caustic model derived in Sec. [3] is used to compute 
the source flux during the caustic exit in order to obtain a syn- 
thetic spectrum at each epoch in a consistent way. In Fig. |2] 
we show the observed and synthetic spectra at two epochs: 9 
July at UT 00h04 and 9 July at UT 22h59 (trailing limb cross- 
ing the caustic). A first analysis reveals good agreement be- 
tween UVES and synthetic spectra for the wings part of Ha 
(6558 - 6561.8 A and 6563.8 - 6567.6 A), whereas a clear 
discrepancy is observed in its core (6561.8 — 6563.8 A). We 
note that the chromosphere is not included in the PHOENIX 
calculations. Al thought it should have a minor influence on the 
broadband limb darkenings, it will have an effect on the core of 
lines like Ha. We therefore divided the analysis in two parts: 
the wings and the core. In order to compare the observed equiv- 
alent widths of the Ha-line with the predictions from our syn- 
thetic spectrum, we apply an overall scale factor to the equiva- 
lent width (1.035 for the wings part and 1.495 for the core part 
of the line), so that the measurements derived from the fiducial 
model match the post-caustic observations. 
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Fig. 2. Upper panel: two UVES spectra (two lower curves) cor- 
responding to 9 July, UT 22h59 (HJD - 2450000 = 2465.47, 
solid line) and 9 July at UT 00h04 (reference spectra, HJD — 
2450000 = 2464.51, dotted line), as well as two computed 
synthetic spectra at the same epochs (the two upper solid 
and dotted curves). Lower panel: fractional difference 5F X = 

2(^04 - F&hsa)/( F obhM + ( lower solid line ) for 

wavelengths in the vicinity of the Ha-line shifted vertically by 

-0.3. Both the observations and the model are shown. On the 

upper right, we show the relative position of the source star at 

each epoch with respect to the fold-caustic shown as dashed. 

The observed and predicted temporal variations of the 
equivalent widths of both the wings and the core of Ha are 
plotted in the upper part of Fig.^ Note at the beginning a small 
short-term decrease in the model, when the leading limb hits 
the caustic, followed by little variation while the source centre 
crosses, but then a marked change when the trailing limb exits 
the caustic. The spectral profile of the core is not well matched 
by the model spectrum as seen in Fig.|2j both the width and the 
depth disagree. Furthermore, the differential variation in flux 
between the centre and the limb is not well reproduced, even 
after the scale factor has been applied. A better fit is obtained 
to the wings. Both parts of the line show larger equivalent width 
variations than predicted by the model. 

Last but not least, in the spectrum corresponding to 10 July, 
UT 02h58 shown in Fig. [3] we clearly note a strong deviation 
in the core of the line but a much smaller one in the wing. For 
the purpose of the following discussion, we define the effective 
photospheric radius of the star as that at which the R band in- 
tensity has fallen to 5 % of the central intensity in R. Using the 
fit done with limb darkening derived from PHOENIX model, 
we infer that the spectrum started when the caustic was at a 
fractional radius p = 0.996 and ended at p = 1.015 while with 
a linear limb-darkening law, the spectrum started at a fractional 
radius p = 1.007 and ended at p = 1.026. 

In both cases, it is clear that the spectrum was taken when 
the caustic was outside the photosphere in R. There is a clear 
emission peak in the core of the Ha absorption line, and this 
can only exist provided there is a temperature inversion in 
the atmosphere (i.e. a chromosphere). The doppler shift of 
the emission core is consistent with material moving outward 
through the giant source star's chromosphere (Ha emission 
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Fig. 3. Same Legend as Fig. |2]but for the pair of spectra 9 
July at UT 00h04. and 10 July, UT 02h58. Notice the strong 
signature of the chromospheric Ha emission. 

line) with a radial velocity of m 7 kms -1 , the signature being 
amplified at the very end of the caustic exit. 
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